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Abstract 

m 

A new phase of temporal evolution of the one-dimensional self- 
gravitating system is numerically discovered. Fractal structure is dy- 
{vq , namically created from non-fractal initial conditions. Implication to 

astrophysics and mathematical physics is discussed. 
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tional many body system. The system represents the motion of parallel 



1 Introduction 



One-dimensional self-gravitating systems is a simplified model of gravita- 



flat mass sheets interacting each other through Newtonian gravity. Since 
its equation of motion can be solved piecewise exactly during each interval 
between collisions, numerical integration of the system is particularly sim- 
ple. Hence the system is used to study long-time behavior and relaxation 
process of gravitational systems pTFo^ , [LSR84| , fMfflj [l'GK97|| . In addi- 



tion, one-particle distribution for canonical and microcanonical ensemble is 
obtained |[Ryb71| . 



Yamashiro et al. studied this model numerically to investigate relaxation 
process of elliptical galaxies [|YGS92|1 . According to the virial theorem the 
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virial ratio V r = 2Eki n /E pot for this model is expected to relax to 1. They 
chose initial conditions with 10~ 3 < V r < 1. Typical feature of the system 
they described is shown in Fig. [I], where a snapshot of distribution in (x, u) 
space (//-space) is shown. We can see a big cluster is formed. 

size-2**10 : init. cond. - water bag: virial ratio-0.1 (step-1200) 
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Figure 1: Particle distribution in /i-space of the model (0) for N = 2 10 ,V r = 
0.1. One big cluster is formed and fractal structure is not present. E=0.25, 
Time=9.375. 




While their main concern was the process of relaxation, they also noted 
that initial conditions with smaller virial ratio ( typically V r < 10 -3 ) did 
not guarantee the relaxation within their CPU time. They also observed 



formation of lumped structure as observed in |[HF67|| , which are quite different 
from the one shown in Fig. |l|. Anomaly in relaxation for small value of V r is 
also pointed out in ||RJM88|| . 

We extend their study to the region V r — > and found that there is a new 
phase of evolution where a fractal spatial structure is spontaneously created. 
We would like to stress that the fractal structure is created even when the 
initial condition is not fractal (typically uniformly random). We also note 
that not all the initial conditions lead to fractal structure. 
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In this letter we report the discovery of this novel phase. In the next 
section we describe the model. We present the numerical results in section 
3, and the final section is for summary and discussions. 



2 Model 

The one-dimensional self-gravitating system (also known as 'mass-sheet model') 
consists of N equivalent flat sheets of constant mass density. The sheets are 
infinitely extended in the y and z direction and are aligned parallel to each 
other and move in the x direction. 

If we suppose the sheets interact with Newtonian gravity, the Hamiltonian 
of the system reads 

N 2 

H = E^ + 27rG ™ 2 EI^-^-l > (!) 

i=l 2m i>J 

where Pi = miij. (ui is the velocity of the sheet i .) We set m = 1/N so that 
the total mass is unity, and 4nG = 1. Time is measured in the unit 



tr 



1/yWjM/L 



where M = mN is total mass, L is the spatial length on which particles are 
distributed initially. t c is called 'crossing time' and is a typical time for a 
particle to traverse the system. 

With this Hamiltonian, the equation of motion reads 

£i = ^ - N l) (2) 

where N R and N L are number of sheets which are on the right and left of the 
sheet i , respectively. 



Eq. (0) indicates that the force acting on each sheet is constant until a 
collision occurs between two sheets. Hence we can exactly integrate the equa- 
tion of motion (^) between every two collisions. This fact greatly enhances 
the usefulness of this model. 

3 Numerical results : structure formation 

Fig. |2] show the evolution of the system (P with iV = 2 15 . Initial condition 
is chosen to be V r = as x% = uniformly random G [0, 1] and «j = 0. 
That is, each sheet is placed randomly with zero velocity dispersion. The 
figures represent particle distribution in (x,u) space (//-space). In the course 
of time evolution we see that fractal structure is formed. In Fig.|^ we show 
box counting dimension of the /i-space distribution shown in the bottom of 
Fig. [2]. Dimension is D = 1.1. 

Fig. ^ is a two-body correlation function £(r) at t — 9.375 in Fig. ^. £(r) 
is defined as 

dP = ndV(l + £(r)) (3) 

where dP is probability to find a sheet in a volume dV at distance r apart 
from a sheet, and n is the average number density. We can see clear power- 
law behavior in £(r) in Fig. 

This power-law correlation is also formed if we modify the initial condition 
as Xi = uniformly random G [0, 1] and Ui = asin27r(a;j — 1/2), as shown in 
Fig. [5[ In this initial condition initial velocity of each particle is uniquely 
determined by the initial position of the particle, hence there is no velocity 
dispersion. 
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"15 : initial condition = uniformly random: velocity fluctuation vamp-0.0 



Figure 2: Formation of fractal structure. Initial condition : Xi : uniformly 
random in [0,1], u { = , N = 2 15 . Time are 2.34375,4.6875,7.03125, and 
9.375. 
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box-counting dimension lor (x,v] space 




Figure 3: Box counting dimension of the //-space distribution shown in the 
bottom of Fig.0. Horizontal axis represents index of division k of the phase 
space, i.e., 2-dimensional /i-space is divided into 2 2k of equal square cells 
whose length i of the side is proportional to 2~ k . Vertical axis represents 
log 2 of number n(k) of the cells which contain at least one particle. Sample 
orbits with the same class of initial condition with different random number 
gives dimension D = 1.1 ± 0.04 for 24 samples. Lines with D = 1 are also 
shown for comparison. 



2-point correlation function : file=n1 5-011 :mesh=10 A 5 2-point correlation function : file=n 15-01 1 :mesh=10 A 5 




Figure 4: two-point correlation function £(r) for t = 9.375 in Fig. |2|. These 
two figures are identical except for the range of vertical axes. Exponent a of 
f oc r~ a is a = 0.20 ± 0.03 for 24 samples. 
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Figure 5: Distribution in (x, w)-space (left) and correlation function £(r) 
(right) for initial condition Xi : uniformly random in [0,1], Ui = 0.1 sin27r(xj — 
1/2) , N = 2 12 . t=6.25. Boxcounting dimension is D = 0.95 ± 0.02 for 24 
samples. Exponent a of the correlation function £ oc r~ a is a = 0.33 ± 0.06 
for 24 samples. 




If we slightly increase the virial ratio from 0, the power-law correla- 
tion is not created. Fig. |6| show the /i-space for an initial condition Xi = 
uniformly random G [0, 1] and w, = uniformly random G [—0.01,0.01]. Here 
we see that fine structure observed in Fig. |2] is no longer observed. Flatness 
of two-point correlation function £ (r) in Fig. [7] supports the impression that 
structure here is not fractal but just made of lumps with constant density. 
This type of state is also observed in | |HF67| and [ |YGS92 |. 



The fractal structure seen in Fig. |2| is not a stationary state. As shown 
in Fig. H the fine structure once created is gradually dissolved into large 
clusters. The reason for the destruction of this structure is not clear now. 

4 Summary and discussions 

In this letter we have shown that there is a new phase in the evolution of one- 
dimensional self-gravitating system where the system creates fractal spatial 
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Figure 6: Example for Vr ^ 0. Fractal structure is not formed. Initial 
condition : Xi : uniformly random in [0,1], Ui = random G [—0.01,0.01] , 
N = 2 15 . Time = 9.375. 
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size=2"1 5 : initial cond. = random(vamp=0.01 ) 
^ n 1 n15-010-mesh-100000 : 1.data" 
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2-point correlation function : file=n15-01 1 :mesh=10 A 5 
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Figure 7: two-point correlation function £(r) for t = 9.375 in Fig. |||. These 
two figures are identical except for the range of vertical axes. 
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Figure 8: Continued from Fig. 0. Time=93.75. 



structure from non-fractal initial conditions. This novel phase appears for 
initial conditions with small virial ratio. The structure is also formed for 
zero velocity dispersion ( i.e., 'cold' condition). Other initial conditions with 
large virial ratio do not lead to this phase. 

As far as we observed, this phase is not a stationary state. Whether or 
not this transiency is a real one or an artifact caused by, for example, our 
computational limitation ( finiteness of number of particles and so on ) is 
not clear now. 

Several studies revealed that two-point correlation function of young stars 



in star-forming regions obey power-law [|Lar95| , [NTHN9q1 . The fractal struc- 
ture we have found may be accounted for an origin of the power-law. Also it 
is a subject of future study if the phase can be found in other systems, such 
as galaxies in the course of evolution. 

Similar structure is observed in a model with expanding universe ||RJF91| , 



TM00| . It will be an interesting subject to examine the effect of expanding 
universe on the structure formation. 
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We would like to stress the importance of the fact that, to create fractal 
structure, it is not necessary for initial conditions to be fractal. Fractal struc- 
ture does not have characteristic spatial scale, nor does Newtonian potential. 
Hence it may not be so surprising that fractal structure can be found in 
systems interacting with Newtonian potential. However the relation between 
these two scale-free objects, Newtonian potential and fractal structure, are 
still unknown. If the initial condition is fractal, since the Newtonian interac- 
tion does not have characteristic length scale, the self-similar property of the 
initial condition would be conserved. In fact, initial conditions with exact 
self similarity (e.g., Cantor set) evolve in time while keeping the self-similar 
spatial structure. On the other hand, we have found that non-fractal initial 
conditions also leads to fractal structure, hence in this case the structure is 
created by the dynamics itself. 

The search for the mechanism why the fractal structure emerges will 
be quite interesting. The emergence is interesting because it is an asymp- 
tote to scale-free behavior. The emergence is also unique in that this is a 
structure formation in conservative system against the tendency to thermal 



relaxation ||KK92| , |KK94 
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cussions, and Yuji Masutomi for useful discussions. We also thank Toshio 
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